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The atmospheric neutrino background for proton decay via p — ► e + iz in ring imaging water Cherenkov 
detectors is studied with an artificial accelerator neutrino beam for the first time. In total, 3.14x 10 5 neutrino 
events corresponding to about 10 megaton-years of atmospheric neutrino interactions were collected by a 1,000 
ton water Cherenkov detector (KT). The KT charged-current single Jl° production data are well reproduced 
by simulation programs of neutrino and secondary hadronic interactions used in the Super-Kamiokande (SK) 
proton decay search. The obtained p —* e + Tt° background rate by the KT data for SK from the atmospheric 
neutrinos whose energies are below 3 GeV is 1.63 ^q 33 ( stat -) -051 ( s y st -) (megaton-year) -1 . This result is 
also relevant to possible future, megaton-scale water Cherenkov detectors. 

PACS numbers: 14.60.Pq,13.15.+g,25.30.Pt,95.55.Vj 



I. INTRODUCTION 

Discovery of nucleon decay would constitute direct evi- 
dence for grand unification of three fundamental forces J3, 
IU, [H and point the way to a new theory beyond the stan- 
dard model of elementary particle physics. Among many 
possible decay modes of nucleons, the decay mode p — > 
e + K° is dominant in a variety of such Grand Unified Theo- 
ries (GUTs) JliH. 

The world's largest (22.5 kiloton fiducial volume) water 
Cherenkov detector experiment, Super-Kamiokande (SK) [7], 
has set a stringent partial lifetime limit of x/B^g+^o > 
5.4 x 10 3 ^__years (90% C.L.) for proton decays into e + and 
Jt° |H, H, [Toll based on the observation of no candidates in 
an integrated exposure of about 0.1 megaton-year (Mtyr). A 
p — ► e + K° signal would be clearly identified in SK as show- 
ering Cherenkov rings corresponding to the positron and two 
gammas from 71° decay, with low net momentum and total 
invariant mass close to the proton mass. The estimated detec- 
tion efficiency for the signal is 40 %, with a background rate 
of about 0.3 events in the integrated exposure predicted by the 
Monte Carlo simulation program (MC). 

Characteristics of the background events from atmospheric 
neutrinos to SK's p — > e + K° search have been studied using 
MC lllll . Charged-current (CC) interactions of atmospheric 
V e with only an electron and single 7t° visible in the final 
state are the dominant source of the background. Neutral- 
current (NC) interactions with only two 7t°s visible in the final 
state are the remaining dominant background. Parent neutrino 
energies between 1 and 3 GeV dominate for both CC and NC 
background events. In this study, SK reports that about 24 % 
of the background events come from neutrinos above 3 GeV. 
The neutrino and final-state nuclear interactions are consid- 
ered to be the dominant uncertainties for the background rate 
estimation. However, the background rate has been estimated 
only with MC without any quantitative systematic error esti- 
mation. 

Some GUT models predict nucleon lifetimes just above the 
current limits, motivating proposals for future, megaton-scale 
SK-type water Cherenkov detector experiments fill [T2I [l3l 
[3,[l5l[i^.[i3l- Since SK and these more massive experiments 
would observe a non-negligible number of background events 
in the future, it is essential to check the predicted background 
rates experimentally. 

In this study, neutrino and secondary hadronic interaction 
MCs used to estimate the background rate in SK proton de- 
cay searches are checked using muon neutrino beam data col- 



lected with a 1,000 ton water Cherenkov detector (KT) in the 
K2K long-baseline neutrino oscillation experiment. The back- 
ground rate for the p — ► e + K° mode in SK and future megaton- 
scale water Cherenkov detectors is determined using the KT 
data. 

While the KT measures muon neutrino reactions rather than 
electron neutrino reactions, the dynamics of pion production 
and re-scattering processes in the oxygen nucleus are identi- 
cal between the two neutrino flavors. Therefore, rare CC V e 
interaction topologies which may mimic proton decay can be 
checked using the KT events with a muon and single 7t° vis- 
ible in the final state. The rate of the atmospheric CC V e back- 
ground can be determined using KT CC single 71° events with 
corrections for differences in neutrino flux and detection effi- 
ciency and the assumption that cross sections and final state 
kinematics for v e and are identical. In addition, the rate of 
the NC background can be determined with two 7t°s visible in 
the final state at the KT. 

The number of Vu interactions in the KT is two orders of 
magnitude larger than previous data sets used for similar back- 
ground studies using other types of detectors JH [H, HHH . 
The K2K beam is well matched to the atmospheric neutrino 
spectrum and samples the energies around a few GeV where 
most backgrounds arise. This is the first result of a proton 
decay background study using accelerator neutrino beam data 
collected with a water Cherenkov detector. 

The outline of this paper is as follows: Section II reviews 
the K2K neutrino beam and the KT detector. Section III-A de- 
scribes the data sample in the KT and validates the neutrino in- 
teraction models. Section III-B calculates the p — ► e + n° back- 
ground rate for proton decay detectors. Finally, Section IV 
concludes and summarizes these results. 



II. EXPERIMENTAL SETUP 

A. K2K neutrino beam 

The wide-band neutrino beam used in the K2K long- 
baseline neutrino oscillation experiment I22I1 was primarily 
muon flavored (about 97.3, 1.3, and 1.5 % for v^, V e , and Vu, 
respectively) with a mean energy of 1.3 GeV. Protons accel- 
erated by the KEK proton synchrotron to a kinetic energy of 
12 GeV were extracted in a single turn to the neutrino beam 
line. The duration of an extraction, or "spill", was 1.1 /jsec. 
The beam was transported to an aluminum target which was 
in the first of a pair of horn magnets. The horn system focused 
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FIG. 1 : The top figure compares the number of total interactions 
in the KT 50 ton fiducial volume for 7.4 x 10 19 protons on target (red 
solid) and the number of total neutrino interactions from atmospheric 
neutrinos in SK (24, 25] for 1 Mtyr exposure (blue dashed for V e +v e 
and blue dotted for V a u+V a u, where all stands for all the neutrino 
flavors). Disappearance of atmospheric v^'s CC interactions due to 
neutrino oscillation is taken into account. The NEUT simulation pro- 
gram [ 26] is used for calculation of total neutrino interaction for each 
neutrino flavor. The bottom figure shows the ratio of neutrino energy 
spectra for the atmospheric neutrino interactions in SK and the K2K 
beam interactions in KT. The meanings of the dashed and dotted 
lines are the same as above. 



positive pions produced in proton-aluminum interactions into 
a 200 meter long decay volume. A beam dump was located at 
the end of the decay volume to absorb all surviving particles 
other than neutrinos. The neutrino beam's profile was mea- 
sured using the distribution of neutrino interaction vertices in 
a muon range detector 112311 . one component of the K2K near 
neutrino detector system. The energy and angular distribu- 
tions of muons from CC neutrino interactions were also con- 
tinuously monitored. 

Figure Q] compares total neutrino interaction (flux x total 
cross section x target volume x time) spectra based on MC 
simulations for atmospheric neutrinos in SK J24l l25ll with to- 
tal interaction spectra for the K2K beam in the KT detector, 
located 300 m downstream of the proton target. These total 
neutrino interaction spectra will be used for the background 
rate calculation shown in Sec. III-B-1. The neutrino spectrum 
measured by the KT covers the same energy range as the por- 
tion of the atmospheric neutrino spectrum which dominates 
the production of proton decay background events. Thus, 
analysis of the KT data allows a controlled study of the neu- 
trino interaction channels and the nuclear re-scattering pro- 
cesses that determine the atmospheric neutrino backgrounds 
to nucleon decay in water Cherenkov detectors. In terms of 
the raw number of neutrino interactions between 1 and 3 GeV, 



the KT v A , data correspond to about 10 Mtyr exposure to the 
atmospheric v e flux at SK. A more complete relationship be- 
tween the exposures will be derived in Sec. III-B-2. 



B. The 1000 ton water Cherenkov detector 

The 1,000 ton ring imaging water Cherenkov detector (KT) 
was located in the K2K near detector hall. The KT was de- 
signed as a smaller replica of SK, using the same neutrino 
target material and instrumentation. The inner volume of KT 
was a cylinder 8.6 m in diameter and 8.6 m in height. This 
volume was viewed by 680 inward-facing 50 cm photomul- 
tiplier tubes (PMTs). The PMTs and their arrangement were 
identical to those of SK, covering 40 % of the detector's in- 
ner surface with active photo-cathode. The software-defined 
fiducial volume for this analysis, based on the reconstructed 
vertex for each neutrino interaction, is 50 tons in a 4 m by 
4 m cylinder at the center of the tank, oriented along the beam 
axis. This fiducial volume is twice the standard fiducial vol- 
ume (25 tons) used in other KT analyses I122T, l27ll in order to 
increase the statistics of the data set. Event selection within 
the fiducial volume results in a pure neutrino sample, with 
negligible contamination from cosmic rays and beam-induced 
muons entering the KT. 

The KT data acquisition system was also similar to that 
of SK. The charge and timing information for each PMT 
hit above a threshold of approximately 1/4 photo-electron 
was digitized by custom electronics modules developed for 
SK 12811 . The KT detector was triggered if there were more 
than 40 hit PMTs in a 200 ns timing window during the spill 
gate. The trigger threshold was roughly equivalent to the sig- 
nal of a 6 MeV electron. The analog sum of all 680 PMTs' 
signals (PMTSUM) was also recorded during every spill by a 
500 MHz FADC, to identify multiple neutrino interactions in 
a single spill. The number of interactions in a spill was de- 
termined by counting peaks in PMTSUM above a threshold 
equivalent to a 100 MeV electron. Only events with a single 
neutrino interaction in the spill, as determined by the PMT- 
SUM peak search, are used in this study. 

The KT and SK share a common library of event recon- 
struction algorithms [29] and detector simulation programs 
based on GEANT3 113011 . Similar detector calibration proce- 
dures are used in the KT [22] and SK, and similar detector 
performances are observed. 



1. Event reconstruction performance 

The KT event reconstruction performance parameters 
which are most relevant to this study are presented in this sub- 
section. These reconstructions are used to select events and 
results of these performance studies are used to estimate sys- 
tematic errors in this study. 

A relativistic charged particle emits Cherenkov light in 
a cone with opening angle cos~'(l/np) ~ 41°, where n ~ 
1.34 is the refractive index for water, aligned with the par- 
ticle's direction, and imaged as a ring by the detector's PMT 
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FIG. 2: Reconstructed - true vertex position of "cosmic ray pipe" 
muons entering at the center of KT. Red crosses (blue boxes) show 
the data (MC). 



grid. Electrons and gammas generate electromagnetic show- 
ers (electrons also suffer multiple scatterings), which produce 
diffused ring patterns, while muons and charged pions create 
sharper ring images, with well-defined edges. 

The interaction vertex for each event is reconstructed using 
PMT timing information. After fixing the vertex, the num- 
ber of Cherenkov rings and their directions are determined 
by a maximum-likelihood procedure. The particle identifica- 
tion (PID) algorithm uses the charge pattern to classify each 
ring as a showering, electron-like particle {e ± , y) or a non- 
showering, muon-like particle (/j*, Tt*). For events with a 
single ring, charge information is used to refine the recon- 
structed vertex based on the PID result. Finally, the particle 
momentum of each ring is determined using charge informa- 
tion, based on the PID classification of the ring. 

Vertex reconstruction in the KT was checked in situ with a 
special device called a "cosmic ray pipe", which is a 25 cm 
diameter PVC pipe with scintillation counters at both ends. 
It was inserted from the top of the tank through a calibration 
access port near the center of the detector's top wall. A coinci- 
dence between the two scintillation counters tagged a cosmic 
ray muon entering the detector at the lower end of the pipe. 
Such a muon has a known direction along the pipe and will 
begin radiating Cherenkov light when it leaves the pipe at a 
known location inside the detector. The "cosmic ray pipe" 
muon data were taken at five vertical positions: 0, ±1, and 
±2 m from the center of the fiducial volume. Figure|2]shows a 
typical reconstructed vertex distribution of "cosmic ray pipe" 
muons. The vertex resolution along the vertical direction was 
10 cm or better, and the central values of the distributions for 
data and MC agree within 4 cm for all test points. The events 
in Fig. |2] have had a vertex refinement based on charge and 
PID information applied to them. The vertex resolution for 



FIG. 3: PID likelihood for cosmic ray muons (shaded histogram) 
and the decay electrons (open histogram). The red (blue) his- 
togram (top (bottom) figure) shows data (MC). Positive (negative) 
likelihood corresponds to muon (electron)-like. 



multi-ring decays without this refinement, such as simulated 
p — > /j + 7I° signal events, is 28 cm. 

The PID performance is checked with cosmic ray muons 
and their associated decay electrons. As shown in Fig. [3] the 
PID likelihood distributions are clearly separated between the 
two types of rings, although note that the momenta of par- 
ticles from the proton decay is typically higher than that of 
Michel electrons. Mis-ID probabilities for 500 MeV/c muons 
and electrons were estimated with Monte Carlo to be 0.5 % 
and 0.8 %, respectively. 

Figure H] shows momentum loss, as estimated by recon- 
structed momentum divided by measured range, for stopping 
cosmic ray muons. The stopping point of a muon is deter- 
mined by the reconstructed vertex of its decay electron. Both 
vertical stopping muons which enter from the top of the tank 
toward the bottom direction and horizontal stopping muons 
entering from one side of the barrel toward the other are shown 
in Fig. |U The largest difference (3 %) between data and MC 
comes from the horizontal muons, and is used in the esti- 
mate of the absolute energy scale uncertainty. The differ- 
ence (1.7%) between vertical and horizontal muons is used 
to estimate the detector asymmetry of the energy scale. The 
directional dependence of reconstructed momentum of Michel 
electrons is also checked and found to be uniform within the 
statistical error (2.4%) of the sample. Figure [5] shows the re- 
constructed invariant mass of NC single-7i° events induced by 
the neutrino beam B27I1 . The events must be fully-contained 
(deposit all of their Cherenkov light inside the inner detec- 
tor) and have two electron-like rings. The K° mass peak is 
clearly observed, but shifted slightly higher than the nominal 
value of 135 MeV/c 2 . This shift is the result of energy deposit 
from de-excitation y's from oxygen nuclei in neutrino interac- 
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energy scale is 1.7%. The absolute energy scale was stable 
within 1 % over the period when data for the present study 
were collected. 



III. DATA ANALYSIS 

All available good quality KT data from the 2000-2004 
K2K running periods were used in this analysis. In this pe- 
riod, 3.14xl0 5 total neutrino events (2.75 xlO 5 events with 
a single neutrino interaction in a spill) in the 50 ton fiducial 
volume, corresponding to 7.4 x 10 19 protons on target (pot), 
were collected. 



A. Validation of neutrino interaction simulations 

Several packages to simulate neutrino interactions have 
been developed OJ. Both the NEUT JH] and NUANCE JH 
simulation programs have been used in SK to estimate the 
background rates for nucleon decay searches. These programs 
have already been used for vario us physi cs analyses and con- 
firmed to reproduce data well ll22il25ll27ll . In this analysis, the 
neutrino-induced /j + k° sample is used to do a careful compar- 



ison of the rates and distributions relevant to the p 
search. 
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data (MC). 



tions as well as a bias in the yy opening angle due to known 
vertex reconstruction bias. Each effect is several MeV and is 
included in MC. The remaining difference between data and 
MC is 4 % by comparing the Jt° mass peaks, and is also used 
to estimate the absolute energy scale uncertainty. Taking these 
measurements together, the absolute energy scale uncertainty 
is determined to be +4 % from the horizontal stopping cosmic 
ray muon and Jt samples, and the detector asymmetry of the 



1. Data and MC event selection criteria 

CC interactions producing only a muon and single Jt° 
visible in the final state ("/J7t " events) can, except for the 
choice of lepton PID, be selected in the KT with nearly the 
same cuts used for the SK p — > e + Tt° search: 

(A) Fully-contained (FC): 

The FC criterion in the KT detector requires that 
no PMT signal greater than 200 photo-electrons 
is recorded, since charged particles leaving the 
detector produce a large signal in the PMT 
around their exit point. 

(B) Two or three identified rings: 

At the neutrino energies relevant to this study, re- 
coiling protons are usually invisible because they 
are below Cherenkov threshold. Events with both 
two and three rings are accepted because, in some 
cases, one of the two gammas from the 7t° decay 
is missed during reconstruction. This can happen 
when the decay is very asymmetric in energy, or 
the rings are too close to resolve. 

(C) PID: 

One identified ring must be muon-like and the 
other one or two rings must be electron-like. 

(D) 71° mass: 

For three-ring events, the reconstructed invariant 
mass of the two electron-like rings must be be- 
tween 85 and 215 MeV/c 2 . 
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The SK p — ► e + K selection criteria are identical, except 
that the outer detector information is used to select the FC 
events, all rings must be electron-like, and events with one or 
more identified muon decays are rejected. 

Finally, the proton decay signal box defines a "p — ► /j7I°" 
sub-sample of "pji°" events consistent with the proton mass 
(about 1 GeV/c 2 ) and Fermi motion in oxygen. The oxygen 
nucleus is modeled with as a relativistic Fermi gas with the 
Fermi surface momentum set to about 220 MeV/c. The selec- 
tion criteria used are: 

(E) Proton decay kinematics: 

Total momentum P tot must be less than 
250 MeV/c and total invariant mass M tot 
must be between 800 and 1050 MeV/c 2 , where: 



• Ptot ■= 

• M tot 

• Etnt 



2 

tot 



2 

-mf, 



• pi is the reconstructed momentum vector of the z'-th 
ring, and 

• mi is the mass of the particle corresponding to the z'-th 
ring inferred from PID (either gamma or muon). 

Figure [6] shows the event display of "p — > pvP" sample in 
the KT real data falling within the proton decay signal box. 
Three Cherenkov rings are clearly visible and well recon- 
structed. 



2. MC data samples 

MC V/j data samples corresponding to 1.1 x 10 20 pot and 
5.4 x 10 19 pot were generated by NEUT and NUANCE, re- 
spectively, and fully reconstructed after full simulation of de- 
tector response using GEANT. These MC samples are used 
to generally study the "pit " events. To extensively investi- 
gate events in the proton decay signal box, additional in- 
teractions were generated with NEUT for 1.4, 5.0, 32.4, 54.0, 
and 35.9 x 10 21 pot in neutrino energy bins 0-2.5, 2.5-3, 3-4, 
4-5, and 5-6 GeV, respectively. This additional sample was 
preselected before reconstruction, using loose cuts on the to- 
tal charge and charge anisotropy from the generated vertex. 
The breakdown of events surviving the selection criteria into 
different reaction channels for the NEUT sample is summa- 
rized in Tbl. U Resonant CC single pion production dom- 
inates both the "pn°" and "p — > ^71°" samples. The quasi- 
elastic contribution arises through production of 7C° by recoil 
nucleon interactions with H2O (also simulated by the GEANT 
and CALOR ll33ll packages). Identification of charged pions 
as muons by the PID is the dominant source of the NC frac- 
tion. 

Figure [7] shows that the rates of events surviving after each 
selection (A-E) for data and both MC samples agree well. The 
simulated data are normalized by the number of total neutrino 
interactions in the 25 ton fiducial volume 112211 in this plot. 




FIG. 6: Event display of "p — > fjltP" sample in the KT real data sat- 
isfying selection criteria (A)-(E). Each hit PMT is represented by a 
small circle whose area is proportional to the PMTs measured charge. 
The small blue triangles show the projection of the reconstructed 
vertex onto the walls. Cyan (red) thick circles show reconstructed 
rings identified by the PID as muon (electron)-like (Two overlapped 
thick circles are identified as electron-like rings). The reconstructed 
M,io, M tot , and P to , for this event are 164 MeV/c 2 , 85 1 MeV/c 2 , and 
226 MeV/c, respectively. 



Event category 



(A)-(D) (%) 



"jUlt sample "p — > //R : sample 



(A)-(E) (%) 

.0.! 



CC quasi elastic 12.2±0.2 (4.3±0.2) 18.6±1.9 

CC single pion from resonance 50.2±0.3 (48.7±0.6) 60.4±2.4 

CCmulti pions 22.2±0.3 (29.9±0.5) 16.1±1.8 

CC deep inelastic scattering 0.6±0.1 (0.5±0.1) 1.0±0.5 

NC 14.8±0.2 (16.5±0.4) 4.0±1.0 



TABLE I: NEUT abundance of different reaction channels 12611 af- 
ter applying event selections (A)-(D) and (A)-(E). Note that hadron 
system invariant mass W is different for "multi pions" and "deep in- 
elastic scattering" modes as 1.3< W <2.0 GeV and W >2.0 GeV, re- 
spectively. "NC" includes all the modes shown for "CC". The errors 
shown are MC statistical. The number in parentheses for (A)-(D) 
shows the population for events with three identified rings only. 



For the data, 24 events remain in the proton decay sig- 
nal box for 7.4 x 10 19 pot, while 21.4 ±3.8 (stat) and 
13.6 ±4.4 (stat) are expected from NEUT and NUANCE, re- 
spectively, for the same number of total neutrino interactions 
in the 25 ton fiducial volume. 

All simulated events in the proton decay signal box were in- 
spected, and about 90 % show a correctly reconstructed, back- 
to-back muon and n° from the CC neutrino interactions. The 
remaining 10 % mostly result from identification of a proton 
as a gamma, or a charged pion as a gamma or a muon. This 
remaining fraction is used to estimate the systematic error on 
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the efficiency difference between SK and KT shown in Tbl. II. 



10 c 



10 = 



at 10" 
c 

0) 

> 

"5 10 3 
a> 

E 

3 2 

c 10 



In the MC, 98 % of "p — > jUJt events found in the signal 
box are from neutrinos interacting with nucleons in the oxy- 
gen nucleus. The remaining 2 % of events are from neutri- 
nos interacting with free protons where a 71° is produced by 
hadronic interaction of the struck nucleon with H2O. 



3. Comparison of kinematics between data and MC 



BCD 

event selection 



FIG. 7: Event rate after application of each selection criterion. The 
event selections (A)-(E) are explained in text. The event rates at D 
and E correspond to "jUJt " and "p — > pK°" samples, respectively. The 
data, NEUT, and NUANCE are shown as red circles, blue squares, 
and blue triangles, respectively. The simulated data are normalized 
by the total number of neutrino interactions in the 25 ton fiducial 
volume. 



Since the background for proton decay searches depends 
not only on the neutrino interaction rate but also on the event 
topologies, it is important to verify that the simulations cor- 
rectly reproduce the kinematics of each particle in the final 
state. Resonant single-pion production, which is the domi- 
nant source of the background, can be characterized by the 
momentum transfer Q 2 and the invariant mass of the hadronic 
system W. Assuming a resonance-mediated three-body reac- 
tion (v^n — > p~A,A — ► pit ) and neglecting the Fermi motion 
and binding energy of the target, the kinematic variables of 
interest can be reconstructed from the observable quantities: 



(E n - p n cos Q % ) {2m„E M - m\ - m 2 ) + (2m n E n - 2E jl E n + 2p jl p- K cos 0^ + m 2 - m 2 ) (m n ~E M + p n cos Q p ) 

(m„ - + pf, cos 9^ - E % + p K cos Q n ) 



and 

Q 2 = 2E V (E M - p^ cos 9^) - m 2 , 

in natural units where p, 9, m, and E are a given particle's 
reconstructed momentum, reconstructed angle from the beam 
direction, mass, and calculated energy, respectively, and 9^ 
as the angle between reconstructed muon and Ji° directions. 
Also, 

_ IntnE,, + W 2 -m 2 n -rrvl 
v 2(m„-E l _ l + p M cosQ l _ l ) 

is the inferred neutrino energy where "tc", "«", and "p" in 
these formulas stand for 7t°, neutron, p~ ', and proton, respec- 
tively. 

Figs.[8]and|9]show the reconstructed hadronic mass W and 
momentum transfer Q 2 of three ring "pn°" events. Uncertain- 
ties in the KT event reconstruction are accounted for in the 
real data points as correlated systematic errors, so the error 
bars shown on each bin are not independent. Although there 
may be a slight excess at small Q 2 in NUANCE, the agreement 
between data and both MC samples is good for both kinematic 
variables within the measurement errors. 



The kinematic variables Q and W do not translate directly 
into the quantities M tot vs. P tot used as the final selection for 
proton decay candidates, so it is important to check the latter 
as well. Figure 10-12 show the distributions of P tot vs. M tot 
for all "pn " events from data, NEUT, and NUANCE, respec- 
tively. 

To simplify comparison between data and MC samples, one 
dimensional projections of M tot for P tot < 250 MeV/c and 
P tot for 800 < Mto, < 1050 MeV/c 2 are shown in Fig.[[3]and 
Fig. Q31 respectively. Note that omitted statistical errors on 
the MC samples are almost the same as those on data. Also, a 
metric quantifying an event's distance from the center of the 
signal box is introduced: L = \JX 2 + Y 2 , where 

X = M tot c 2 - 938 MeV 




P tot c - 200 MeV (P tot > 200 MeV/c) 
(P M < 200 MeV/c) 



The L distributions for data and MC samples are plotted in 
Fig-OS The data and both simulations agree well, within mea- 
surement errors, for all values of M tot (Fig.[T3~l>. P tot (Fig.[T4l>. 
and L (Fig. IT5l l. including the tails of distributions where the 
proton decay signal box is located. 
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reconstructed hadronic invariant mass (MeV/c 2 ) total invariant mass (MeV/c 2 ) 



FIG. 8: Reconstructed hadronic invariant mass W for three ring 
"jUJt " events, assuming CC single 71° production via resonance. Red 
crosses show the data with statistical and total measurement errors. 
The relatively large systematic errors in this plot are correlated be- 
tween bins, and arise from the absolute energy scale uncertainty. The 
NEUT (NUANCE) prediction is shown by the solid (dashed) blue 
histogram. The hatched histogram shows the distribution of NEUT 
events that originate from resonant CC single pion production. Both 
MC samples are normalized to the data by number of entries. 




reconstructed Q 2 (GeV 2 ) 

FIG. 9: Reconstructed momentum transfer Q 2 for three ring "jUir " 
events, assuming CC single 7t production via resonance. Red 
crosses show the data with statistical and total measurement errors. 
The NEUT (NUANCE) prediction is shown by the solid (dashed) 
blue histogram. The hatched histogram shows the distribution of 
NEUT events that originate from resonant CC single pion produc- 
tion. Both MC samples are normalized to the data by number of 
entries. 



FIG. 10: Total momentum P tot vs. total invariant mass M, or for "/jit 
events from the KT data (7.4 x 10 19 pot), with the proton decay 
signal box superimposed. 




'o 200 400 600 800 1000 1200 

total invariant mass (MeV/c 2 ) 



FIG. 11: Total momentum P tot vs. total invariant mass M, or for "/jit 
events from NEUT (1.1 x 10 20 pot), with the proton decay signal box 
superimposed. 



In summary, these comparisons show good agreement be- 
tween data and MC samples; hence, the modeling of neutrino 
and final-state nuclear interactions relevant to p — > e + 7i° back- 
grounds appears to be well supported by the KT data. 
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FIG. 12: Total momentum P to t vs. total invariant mass M, or for "/m 
events from NUANCE (5.4 x 10 19 pot), with the proton decay signal 
box superimposed. 



B. Determination of the background rate to the p — * e + 7t 
search 



FIG. 13: Total invariant mass M, or for "pTl events with total mo- 
mentum Ptot < 250 MeV/c. Red crosses show the data with sta- 
tistical and total measurement errors. The NEUT (NUANCE) pre- 
dictions are shown by a solid (dashed) blue histogram. Both MC 
samples are normalized to the data by number of entries. The green 
arrow shows the range accepted by the proton decay signal box. 



1. General method 



Predicting the background rate for the p — > e + n search re- 
quires extrapolation from the KT v A , beam data to all flavors 
of the atmospheric neutrino flux, and re-weighting the K2K 
neutrino spectrum to match the atmospheric neutrino flux and 
spectrum. The CC background event rate from v e is obtained 
using lepton universality, the assumption that the cross sec- 
tions and final state kinematics forVg and are nearly iden- 
tical. The NC background rate is directly measured. The 
anti-neutrino background rate for all flavors is estimated by 
rescaling the data to the predicted total event rates for each 
anti-neutrino. Anti-neutrino final-states tend to have smaller 
Q 2 due to the different fraction of neutrino energy transferred 
to the hadronic system in inelastic reactions and therefore tend 
to have a larger momentum imbalance between the lepton and 
hadronic systems. Simulated data samples confirm that anti- 
neutrinos are less likely than neutrinos to fall within the signal 
box, so the rescaling procedure is conservative. 

The expected event rate of the p — > e + n° background for 
general water Cherenkov detectors, N, can be expressed as: 



10' 



jj 10 

E 

3 
C 



^ -t-i-t'--t-H - 



100 200 300 400 500 600 700 800 900 1000 

total momentum (MeV/c) 



N = n cc ■ R% c ■ RC C 



„NC pNC nNC 



(1) 



where the 1st (2nd) term corresponds to the background rate 
coming from the CC (NC) neutrino interactions, n is the ob- 
served events in the proton decay signal box at the KT, and 
is a ratio of total neutrino interactions (flux x total cross sec- 
tion x target volume x time) between the atmospheric neutri- 
nos at the proton decay detectors and the K2K muon neutri- 



FIG. 14: Total momentum P tr 



for 

„2 



"it/it"" events with total invariant 
mass 800 < M tot < 1050 MeV/c z . Red crosses show the data with 
statistical and total measurement errors. The NEUT (NUANCE) pre- 
dictions are shown by a solid (dashed) blue histogram. Both MC 
samples are normalized to the data by number of entries. The green 
arrow shows the range accepted by the proton decay signal box. 
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600 800 

L (MeV) 

FIG. 15: L (as defined in the text) for "fm°" events. Red crosses 
show the data with statistical and total measurement errors. The 
NEUT (NUANCE) predictions are shown by a solid (dashed) blue 
histogram. Both MC samples are normalized to the data by num- 
ber of entries. The green arrow indicates the range accepted by the 
proton decay signal box. 



nos at the KT. In using Eqn.Q] V e +V e and V„//+V„// (where all 
stands for all the neutrino flavors) are used for the atmospheric 
neutrino flux for the CC and NC background measurements, 
respectively. R§ is essentially an energy-dependent flux cor- 
rection factor, multiplied by an overall scale factor based on 
the respective target masses and exposure. Finally, R £ is the 
ratio of detection probabilities for the background events at 
the proton decay detectors and the KT. Any difference in de- 
tection efficiency between the "p — > e7t°" sample in the pro- 
ton decay detectors and the "p — > pK°" sample in the KT is 
reflected in R^ c . 

Since the numbers of "p — > p%°" events from NEUT and 
NUANCE at KT are statistically consistent with each other, 
one (NEUT) of the models is chosen for the background rate 
determination in this study. 

As shown in Sec. III-A-2, the remaining number of "p%°" 
events in the proton decay signal box, n , is 24 for 7.4 x 
10 19 pot. 

The number of NC background events, the NC interactions 
with only two 7i°s visible in the final state, is measured with 
the KT by applying the same event selection criteria used in 
the SK p — > e + n° search. The events must be fully-contained, 
and have two or three rings, all of which must be electron-like. 
For three-ring events, there must be a reconstructed 7t° mass 
between 85 and 215 MeV/c 2 . The total mass of the event must 
satisfy the conditions 800 < M tot < 1050 MeV/c 2 , and the P tot 
< 250 MeV/c. 

In the SK data analysis, a decay-electron cut is also applied 
to the data and MC samples. The reconstruction software used 
at the KT during K2K running did not implement a decay- 



electron cut. However, it was found that the NC background 
rate was so small after the previous cuts that there were no 
data candidates left to check for muon-decays. For the KT 
MC "p — ► en " sample, a cut was applied using the MC truth 
information assuming the same decay-electron finding effi- 
ciency as SK. With the exceptions of the decay-electron cut 
and method for FC event selection, which in SK employs the 
outer detector, the selection criteria applied for these two de- 
tectors are identical. 

No candidates of the "p — > e7t°" events are found in the KT 
data (n NC = for 7.4 x 10 19 pot), while the expected number 
of NC events is 0.30±0. 13 (stat.) for the same number of total 
neutrino interactions in the 25 ton fiducial volume. All the 
MC "p — > e7t°" events from the NC interactions are found to 
have only two 7i°s visible in the final state. These events were 
NC single pion resonance events with an extra pion produced 
by hadronic interactions in the water. 



2. Application to Super-Kamiokande-I 

In order to calculate the background to the proton decay 
search at the SK-I experiment [9, 10] the interaction rates are 
first re-weighted to the atmospheric flux in SK ll24l l25ll and 
the SK fiducial volume of 22.5 kiloton and used to calculate: 



R CC(NC) = 



^energy bins CC(NC) CC(NC) 



energy bins ^CC{NC) 
j 



Hj m 



(2) 



where m cc ^ NC ^ is the expected number of "p — > p(e)%°" events 
from CC(NC) interactions in the KT for 7.4 x 10 19 pot and 

CC(NC) 

n, is the expected ratio of V e +V e (v„//+V„//) interactions 



in SK for 1 Mtyr to the interactions in the KT. r, 



CC(NC) 



and 



m CC{NC) are s }j 0wn m pjg [JJ (bottom) and Fig. Q2J respec- 
tively. 

The background rate is estimated only for E v < 3 GeV in 
this study since the expected number of "p — > p(e)n°" events 
in the KT is small (<10 _1 ) above 3 GeV as shown in Fig. [16] 
The obtained results for E v < 3 GeV are R^ c = (15.9 MtyrT 1 

and R% c = (4.5 Mtyr) -1 . For example, R^ c = (15.9 Mtyr)" 1 
means that KT Vu data corresponds to 15.9 Mtyr exposure of 
the atmospheric neutrino data of the CC background events at 
SK. The difference between R^ c and R^ c mostly comes from 

the difference between r^ c and r^ c . The fraction of back- 
ground events above 3 GeV in a search for p — * e + n° at SK is 
reported to be about 24 % using MC full . 

The efficiency ratio for the background events, R E , 
can be estimated by using proton decay signal MC (note 
that the majority of the background events really have visi- 
ble Cherenkov rings and their kinematics are consistent with 

proton decay signals as shown in Sec. III-A-2): R £ = e 
I £ CC(NC)^ ^ere e = 0.40 is the efficiency for p — > e + %° MC 
events at SK and e cc ( NC ) = 0.37(0.34) are the efficiencies for 
p — >/j + 7l° and p — > e + n° MC events in the KT. The largest 
source of inefficiency for both detectors is due to final-state 
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neutrino energy (GeV) 



FIG. 16: Predicted number of CC (m cc (i) shown in top) and NC 
(m NC (i) shown in bottom) "p —* ,u(e)n " events from NEUT in the 
KT for 7.4 x 10 19 pot, in bins of true parent neutrino energy i. The 
errors shown are MC statistical. 



Error source 


error on rf c (%) 


KT vertex reconstruction 


3 


KT absolute energy scale 


3 


KT detector asymmetry of energy scale 


3 


KT FC event selection 


1 


KT ring counting 


9 


KT ring direction 


3 


KTPID 


9 


KT NC fraction in "p — > /J3t°" events 


+0/-4 


Sub total 


+ 14/-15 




Error source 


error on R^ C{NC) (%) 


K2K flux shape 


3(6) 


K2K beam v e contamination 


<K5) 


KT MC normalization 


4(4) 


KT MC statistics 


5(47) 


atmospheric v flux [24, 25] 


10(10) 


total v cross section ratio 


15(15) 


lepton universality 


<K0) 


atmospheric anti-v rate 


+0/-16(6) 


Sub total 


+257-29(54) 



pion interactions in the oxygen nucleus. The efficiencies for 
free protons are 0.80 and 0.77 for p — > e + n° MC events at 
SK and p -^p + n° MC events in the KT, respectively. The 
15 % (=0.34/0.40) difference between SK and KT efficiencies 
for p — > e + 7t° is mainly due to the different FC selection cri- 
terion. 

Finally, using the numbers above the atmospheric neutrino 
background rate for the SK p 
to be: 



Error source 



error on R £ (%) 



e + n° search, N, is determined 



N {Mtyr' 1 ) 



{ 1.63 ± 0.33 (stat.) t° fi {syst . ) } (CC) 
-{0.00 + 0.26{stat.)+0.13(syst.)}(NC) 
= lM±™l(stat.)t f 1 (syst.). 



Table [TT] summarizes the systematic errors on each param- 
eter in Eqn. (fTJ. Uncertainties from the performance of the 
KT reconstruction are reflected in the error on n . The error 
is estimated by evaluating the effect of reconstruction uncer- 
tainties on the number of p — ► p + n° MC events in the proton 
decay signal box. To account for the NC contamination in the 
"p — > pn°" CC events (Tbl. U), an asymmetric error of (+0/- 
4%) is assigned to n cc . It was not necessary to estimate 
the errors on n NC as there was no observation of "p — > e7t°" 
events in the KT. The expected event rate in this case is much 
smaller than that from the CC events, which dominates the 
background. 

One (15 %) of the dominant errors on R^ c comes from un- 
certainty of total neutrino interactions with atmospheric neu- 
trino flux at SK. The cross section systematic errors consid- 
ered are 10% on the axial vector masses for quasi-elastic 
and single -pion production, 10% on the total cross sections 
for quasi-elastic and single-pion production, 5 % on the total 



Efficiency difference between SK and the KT 



11(7) 



TABLE II: Summary of systematic errors for measurement of the 
background rate to the p — > e + n° search. 



cross section for deep-inelastic scattering, with and without 
correction of the nuclear structure functions for deep-inelastic 
scattering (Bodek & Yang correction Il34ll ). and 20 % on the 
ratio of NC/CC cross sections. The other dominant error (+0/- 
16 %) on R^ c comes from differences in kinematics of the 
final-state particles between neutrinos and anti-neutrinos. The 
error on R^ c (R^ ) is estimated by comparing the number of 

simulated CC single-7t° (NC two-7t°) events near the signal 
box. The error on R^ c from the lepton universality assumption 

is estimated by comparing the number of CC single-7t° events 
near the signal box using simulated Vu and v e interactions and 
is found to be negligibly small. The largest error (47 %) on 
R^ c arises from the statistics of simulated "p — > eK°" events 

(m NC in Eqn. ©) in the KT. The KT MC data are normal- 
ized by the number of total neutrino interactions observed in 
the 25 ton fiducial volume. The estimate of the error on the 
fiducial volume is 4% 12211 . 

CC( NC) 

The systematic error on R £ is estimated by taking into 
account the differences of PID performance between SK and 
the KT and the mis-identification probability of charged pions 
or protons as either muons or gammas in the KT MC sample. 
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3. Discussions 

A total of 2.14 (=1.63/0.76, where 0.76 is fraction of the 
background events below 3 GeV at SK lll IIP atmospheric neu- 
trino background events are expected for 1 Mtyr for SK-type 
water Cherenkov detectors. If the proton lifetime x/B p ^ e+1 p 
were 6x 10 34 years, the proton decay signal rate would be al- 
most equivalent to the rate of the background events. There- 
fore, it will be important to reduce the number of the back- 
ground events in order to make a clean discovery of pro- 
ton decay in future experiments. If necessary, the remain- 
ing background can be further reduced to improve the sig- 
nal/background ratio. For example, Fig. [14] suggests that the 
background rate could be reduced by an order of magnitude by 
applying a tighter momentum cut (P tot <100 MeV/c), while 
the signal efficiency would be reduced only by a factor of 2.3 

With the present SK exposure (0.092 Mtyr), the expected 
background rate for E v < 3 GeV is 0. 15 t° Q i (stat) +°;° 4 (syst) 
events. This result is consistent with the observed lack of can- 
didates in the SK experiment ifoL fToh . 

According to Tbl.[II] the background rate determination can 
be improved by reducing uncertainties of the atmospheric neu- 
trino flux, and uncertainties of total and anti-neutrino interac- 
tion rates with the atmospheric neutrino flux in future. 

Note that the background rate for the p — > /j + k° mode can 
be determined in the same way just by replacing the atmo- 
spheric V e +V e flux with V^+V^ flux for the CC background and 
the detection probability of the "p — > en " sample with that of 
the "p — > /J7t°" sample in the proton decay detectors. 

IV. CONCLUSION 

The atmospheric neutrino background to searches for p — > 
e + K° has been experimentally studied using an accelerator 
neutrino beam and a water Cherenkov detector for the first 
time. The K2K KT detector, with the same target material 
and detection technique as SK, accumulated data equivalent 
to atmospheric neutrino exposures of 15.9 and 4.5 Mtyr for 
the CC and NC background events, respectively. 

Two neutrino interaction simulation programs, NEUT and 
NUANCE, were evaluated and found to reproduce accelerator 
neutrino beam interactions in water within the measurement 



errors. Both should therefore reliably model the atmospheric 
neutrino background to p — ► e + n° searches. 

This is the first determination of the background rate us- 
ing data for the SK-type water Cherenkov detectors. We mea- 
sured the rate of neutrino and final-state nuclear interactions of 
the background events. To measure the background rate from 
the KT data themselves, almost identical selection criteria as 
for p — ► e + 7t° searches were applied. A total of 24 "proton 
decay"-like events were identified. A re-weighting procedure 
was applied using the atmospheric neutrino flux and detection 
efficiency (40 %) for the SK proton decay search, and the re- 
sulting background rate for E v < 3 GeV (where about 76 % of 
the background events is expected in SK according to the MC 
study) is 1.63+[J;^(stat)+JJ;^ (syst) Mtyr" 

This experimentally determined background rate is consis- 
tent with no candidates being observed over the 0.1 Mtyr ex- 
posure of the SK-I experiment. 

This result shows that about two background events per 
year would be expected in possible future one-megaton-scale 
detectors. Assuming a finite proton lifetime by an order be- 
yond the present limit, the rate of signal would be similar to 
the expected background rate, both in these proposed detec- 
tors as well as in a still-running Super-Kamiokande. There- 
fore, further reduction of the background events will be crucial 
in the future. 
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